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Zero-temperature responses of a 3D spin glass in a field 
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We probe the energy landscape of the 3D Edwards- Anderson spin glass in a magnetic field to test 
for a spin glass ordering. We find that the spin glass susceptibility is anomalously large on the lattice 
sizes we can reach. Our data suggest that a transition from the spin glass to the paramagnetic phase 
takes place at B c ~ 0.65, though the possibility B c — cannot be excluded. We also discuss the 
question of the nature of the putative frozen phase. 
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Ising spin glasses |0 have been studied intensively for 
over two decades; nevertheless, even the basic issue of 
the nature of the phase diagram is still unsettled. There 
are two main schools of thought in the long standing de- 
bate over whether a spin glass phase can exist in the 
presence of a magnetic field B. In the mean field pic- 
ture B the spin glass phase exists up to some critical 
field value Bat{T), i.e., up to the so called "de Almeida- 
Thouless" line H (AT) where replica symmetry breaking 
(RSB) arises. In the droplet picture QJ^], any non-zero 
magnetic field kills the spin glass ordering and makes the 
system paramagnetic. 

Experimental evidence in favor of each school of 
thought has been claimed ||J7|], but no consensus has 
emerged. On the computational side, the approaches us- 
ing Monte Carlo have been hindered by the large finite 
size effects present in d — 3 [p|-p0[. Because of this, 
dynamical out-of-equilibrium approaches have been used 

{e.g., 



nM for d = 4 and 1 13 LJ for d 



3); those works 
hint at the existence of a phase transition in field, but 
it is preferable to have a direct test using a system in 
equilibrium. For such a test, one should stay away from 
the critical point at B = 0. To achieve this, we focus in 
this work on the T = line in the B-T phase diagram; 
then the spin glass ordering can be tested through the 
behavior of excitations above the ground state. 

We begin by discussing possible scenarios for the spin 
glass ordering when B ^0. In our first approach, we ex- 
tract low energy excitations and investigate their size as 
a function of L and of B; we also consider some of their 
topological properties. In a second approach, we compare 
ground states with periodic and anti-periodic boundary 
conditions, and extract a "magnetic penetration length'' . 
We use extrapolations and finite size scaling to estimate 
the critical field B c where the system becomes paramag- 
netic. Our results suggest B c rs 0.65. 

The model and spin glass ordering - We consider 
the EA Hamiltonian on a 3D L? periodic cubic lattice: 



HjdSiV^-^JijSiSj-B^Si 



(1) 



<i]> 



In this Edwards-Anderson (EA) model, the sum is over 
all nearest neighbor spin pairs, Si — ±1, and B is the 
strength of the magnetic field. The quenched couplings 
Jij are independent random variables taken from a Gaus- 
sian distribution of zero mean and unit variance. At 
B = 0, it is widely believed that the low T behavior 
of this system is characterized by a frozen but random 
ordering of the local magnetizations: (Si) ^ 0, where (.) 
is the thermal average. When B > 0, there is no up-down 
symmetry and the fact that (Si) / is not relevant. In- 
stead, one relies on the correlation length £: £ = 00 in the 
spin glass phase, whereas £ is finite in the paramagnetic 
phase. 

We study the system at T — and ask what signature 
corresponds to a £ = 00 behavior of T > system. Given 
the ground state in the presence of B, consider droplet 
excitations, i.e., excitations that arc of minimal energy 
at fixed "size" (taken generally as the "radius" r of the 
connected cluster of spins that are flipped) . Let be the 
exponent B giving the characteristic excitation energy of 
these droplets: E rs r e . If the distribution of the rescaled 
energies, Pd(E/r e ), has a non-zero density at zero argu- 
ment, then there are enough thermally activated droplets 
of all sizes to make the correlation length infinite for low 
enough temperatures. We have three main scenarios: 

(1) The "mean field" scenario: if B stays below some 
critical value Bat there is RSB ; for B > Bat the be- 
havior is paramagnetic. Because of RSB, spin-spin cor- 
relation functions do not cluster when B < Bat'- any 
reasonable definition of £ leads to £ = 00 below Bat- 

(2) The "scaling/droplet" scenario: as soon as B > 
the system is paramagnetic. Motivation for this scenario 
comes from the large field limit: there, 8 = 3 because ex- 
citation energies grow linearly with the number of spins 
flipped. Based on fixed-point arguments, one then ex- 
pects = 3 and Pd(0) = for all B > 0. Because of this 
last property, large scale droplets cannot be thermally 
activated in the presence of a field (or equivalently when 
there is an extensive magnetization). 

(3) An "intermediate" scenario: the spin glass order 
survives at positive B, but there is no RSB and thus no 
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FIG. 1. Mean volume of the droplets extracted 
one-spin flip method as a function of L. From top 
torn: B - 0.1, 0.2, ■ • •. In the insert B = 0.8, 1.0 and 
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AT line. 9(B) may be arbitrary while B c is denned as 
the field where P<z(0) goes to 0, separating the regime of 
£ = oo from the one where £ < oo. Since this scenario 
has large scale thermally activated droplets co- existing 
with an extensive magnetization, we will refer to it as 
the MAD (Magnetization And Droplets) scenario. 

Droplet sizes - In our first approach we generate 
low energy clusters of spins according to a one spin flip 
method. Given the ground state Co of H, we randomly 
choose a spin Si and force its orientation to be opposite 
from what it is in the ground state. Then we recom- 
pute the new ground state C given this constraint. (The 
numerical algorithm for computing the ground states is 
described in |lq] .) The difference between C and Co is a 
connected cluster of spins whose volume we denote by V] 
it is the lowest energy cluster among all those containing 
Si . The mean cluster size (V) is analogous to the spin 
glass susceptibility xsg — L 3 ((q — (q)) 2 ), where q is the 
standard spin overlap. A diverging (V) as L — > oo corre- 
sponds to a spin glass phase, while a bounded (V) is the 
signature of a paramagnetic phase. 

For each value of B (0.1 < B < 1.2, in steps of 0.1) we 
have generated 3000 disorder instances for L = 6,8, 10, 
and 1000 for L = 12. In figure pi we display the mean 
volume (V") of the droplets generated by this one-spin 
flip method. (The spin Si selected for flipping is chosen 
at random amongst the 25% spins with the largest local 
fields; this enhances the signal but does not affect our 
conclusions.) At B > 1 we see signs of (V) saturating 
when L grows: these large B values are in the paramag- 
netic phase. For smaller values of B, (V) grows signifi- 
cantly with increasing L. In the droplet picture there is 
a magnetic length £b which acts as a cut-off: (V) will 
grow as a power of L until L ~ £g, and thereafter it will 
saturate. Furthermore, this cut-off length grows as B de- 
creases, Ib ~ (Tf/B) - 2— . On the contrary, in the MAD 
and mean field scenarios (V) diverges for all B < B c . 

Consider for example the data at B = 1.0. If we use 
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FIG. 2. Fraction of sponge-like excitations versus B. 



the estimates on the largest lattices 1 16 available to date 
(T w 1.78 and w 0.19), we have £ B {B = 1.0) w 2.13. 
(V) should be roughly the cube of this number, while we 
find (V) — 52 at L = 12. Such a large value of (V) is 
unexpected in the droplet picture. Unfortunately, since 
the values of T and 9 have large uncertainties, one can- 
not be more quantitative than that. In fact, the curves 
displayed in figure [l] could be interpreted either by saying 
that for B < 0.8, Ib is comparable to our largest lattice 
sizes, or that (V) diverges as L — > oo even when B is not 
too small, for instance for B < 0.6. To further test the 
different scenarios, we follow |17| and focus now on the 
topology of our excitations. 

Fraction of sponge-like excitations - In the mean 
field (RSB) picture, macroscopically distinct valleys dif- 
fering by energies of 0(1) arise with positive probability 
when B < Bat and with zero probability for B > Bat- 
(When the volume is finite, this last probability should 
go to zero exponentially in L.) These valleys differ by 
excitations that span the whole system and that wind 
around the lattice. We say that the cluster associated 
with one of our excitations is sponge-like if both it and 
its complement wind around all three directions (x, y, z) 
of the lattice. This motivates our measuring the fraction 
f(B, L) of the excitations generated by our one-spin flip 
method that are sponge-like. (For low field values, we 
observe that the flipped cluster sometimes contains more 
than 0.5L 3 lattice sites; since such events disappear as 
L — > oo in all three scenarios, we have excluded them 
from our analysis in order to reduce the finite size correc- 
tions.) Let f*(B) be the large L limit of f(B, L). f*(B) 
is an order parameter for RSB: it is zero for B > Bat 
and positive for B < Bat- 

In figure g we show our values for the fractions f(B, L). 
As expected at large B, f goes to zero with increasing 
L. If there is RSB, the large L curves will converge to 
a limiting function that intersects the x axis at Bat > 
0. On the contrary, in the other pictures where replica 
symmetry is not broken, the limiting curve is f*(B) = 



for all B > 0. 

What can be said from these data within the RSB sce- 
nario? The order parameter should behave as f*(B) ss 
(Bat — B) 13 with j3 less or equal to its mean field value 
of 1. Thus we expect /* to be convex for B not too 
close to 0. Now if we make the reasonable assumption 
that f(B,L) converges to /* monotonically, then we get 
a bound on Bat by taking the tangent to the L = 12 
curve; this leads to Bat < 0.65. Of course, we can be 
more realistic and say that the curves for different L drift 
towards /* in a smooth way. If we do that and extrapo- 
late the curves by eye, we find Bat < 0.4. Obtaining a 
less subjective estimate requires using finite size scaling 
and parametrizing /*, but our small range in L makes 
such an attempt futile. 

Consider now the data from the point of view of the 
two scenarios with no RSB. We have performed the fi- 
nite size scalings f(B, L) = L"0A , F[(B C -B)L 1 /' J ] where 
B c = in the scaling/droplet scenario and B c > in the 
MAD scenario. For each putative value of B c , we first 
adjust fl/v so that F is nearly L independent at B = B c , 
and then we adjust v for the curves to superpose as well 
as possible. For both B c > and B c = the data col- 
lapse is far from perfect, so we cannot use the chi squared 
for a quantitative analysis. However, the quality of the 
superposition can be judged by eye. For B c w 0, the 
superposition is reasonably good and leads to v = 1.25 
and j3 ~ (recall that sponges seem to arise with a fi- 
nite probability in zero field). As B c is increased, the 
superposition first becomes less good but then becomes 
better again at intermediate values with a local optimum 
at B c w 0.65; there we find v w 0.8 and (3 w 2.0. With- 
out a reliable parameterization of the finite size effects 
and larger lattice sizes, it is not possible to go beyond 
these qualitative conclusions. 

Periodic- anti-periodic computation - In our second 
approach we perturb the couplings Jij and probe how this 
affects the ground state. First we compute the ground 
state with periodic boundary conditions. Then we change 
the sign of all the Js cutting a given vertical plane II; 
this amounts to applying anti-periodic boundary condi- 
tions across this plane. Finally we re-compute the ground 
state for this new system. In contrast to the zero field 
case, gauge invariance is broken here and the plane II 
where the Js are reversed is measurable. Our reason for 
choosing this particular perturbation is that it is trans- 
lationally invariant in two directions; this leads to good 
statistics on the observables. 

Given this "perturbation" , how does the difference be- 
tween the two ground states look? At high fields, a few 
spins will be flipped in the immediate neighborhood of 
the plane II. As B decreases, the region affected by the 
perturbation will broaden; one can introduce a "mag- 
netic penetration length" £p as a measure of the region's 
width. In the paramagnetic phase, £p is finite, while it 
is plausible that it grows with the system size L in a 
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FIG. 3. Inset: the raw data for £p versus the field B. Main 

plot: £p rescaled by L 0,42 versus rescaled distance from B c . 



spin glass phase. It is possible that in a droplet picture 
£p = £b , but we do not develop this issue here. 

To define £p from the perturbed region, we measure 
the spin-spin overlap of the two ground states and then 
average in each plane parallel to II. Let Q(d) be this 
overlap, with d the distance from II. In the paramagnetic 
phase Q(d) will go to 1 at large d, and the approach 
to this asymptote should go as exp(— d/£p). If on the 
contrary £p = oo, the asymptote may be different from 
1, and the approach to that value can be a power law in 
d. For each L and B we extract £b by fitting Q(d) to 
the form 1 + Aexp(— d/^s). Above B c , this should lead 
to the measurement of the true £b, while below B c the 
extracted £p will be an effective length that will grow 
with L. (In these measurements, we used 40000 disorder 
samples at L = 6 and 20000 at L = 8 for each B; for 
L = 10, we used 5000 samples at the smallest B value 
and up to 12800 on the largest one, while we used from 
1000 to 7000 samples for L = 12.) The fits performed 
using all the distances gave results very similar to the 
ones using only the four (symmetrized) data points at 
the largest values of d. 

We show the results of our best fits in the inset of 
figure pi £p saturates at high fields, while it increases 
with L at low field. We also show in the figure a finite 
size scaling plot in which we have set B c = 0.65; £p has 
been rescaled by L~ QA2 and (B - B c ) by L 025 . The data 
collapse is good, but similarly good results are obtained 
for all values of B c < 1.0; the problem is that when data 
do not have a big dynamic range, finite size scaling can 
almost always be made to work. Thus for this method 
we have no precision on the estimate of B c . 

We have also used other methods to estimate £p from 
the set of spins that are flipped in the periodic-anti- 
periodic transformation; these include the set's mean dis- 
tance to II, the mean extension of its interface, etc... All 
these measures lead to £p comparable to L when B < 0.7 
for our range of Ls. It is thus difficult for us to say 



gc 

1 




FIG. 4. Fraction of sponge-like clusters in the periodic 
anti-periodic approach. The insert gives the best finite size 
scaling fit (B c = 0.62, v = 1.0). 



study leads us to suggest that Bat < 0.4 if it is in- 
deed positive. This very low value is compatible with 
the findings of [[13]]. However, confirming such a value 
via equilibrium measurements will require larger lattice 
sizes than we can handle at present. 
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whether £p diverges when B reaches B c > or only when 
B — > 0, but this is not unexpected, and has been a serious 
shortcoming in all attempts to detect a spin glass order- 
ing in a field. Because of this, we now again appeal to 
topology: we measure the proportion of instances where 
the flipped spins form a sponge-like cluster. 

The raw data for the fraction of sponges are plotted 
in figure 0. We observe a crossing point, suggesting a 
first order jump in this fraction at B c rs 0.65. This is 
probably our clearest evidence for the existence of a spin 
glass phase with B c > 0. To perform finite size scaling, 
we take the fraction of sponge-like events to be a function 
of [B — B c )L l / u ] the best data collapse is obtained for 
B c s; 0.62 and v w 1.0 and is displayed in the figure's 
insert. (Note that the data at low B goes to the limiting 
curve only slowly when L increases.) These data are 
compatible with a B c w 0.65; however, since the crossing 
points drift a bit, a smaller value for B c cannot be ruled 
out. 

Discussion - We have introduced and exploited a 
new approach for testing whether a spin glass ordering 
arises in the presence of a magnetic field. The major 
advantage of our method compared to Monte Carlo is 
that since we work at T = 0, we are far from the B = 
critical point, T c w 0.95. Our main conclusion is that 
very plausibly in three dimensions, a spin glass phase 
survives up to a critical field B c sw 0.65; beyond that, the 
system becomes paramagnetic. 

Our proposed value is several times smaller than the 
mean field value p8[ , Bj£f « 2.1. Such a low value for 
B c makes it difficult to establish without ambiguity that 
B c > 0. This problem is analogous to the difficulty of 
showing that T c > in d = 3 spin glasses; to give plau- 
sible evidence, finite size effects must be under excellent 
control. Clearly, that is not yet the case here, but the fact 
that our different approaches lead to consistent values of 
B c gives some strength to our findings. 

Finally, there is the question of a possible replica sym- 
metry breaking transition at Bat > 0. Our numerical 
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